Magnon BEC and various phases of 3D quantum helimagnets under high magnetic field 
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We study high-field phase diagram and low-energy excitations of three-dimensional quantum helimagnets. 
Slightly below the saturation field, the emergence of magnetic order may be viewed as Bose-Einstein conden- 
sation (BEC) of magnons. The method of dilute Bose gas enables a quantitative analysis of quantum effects 
in these helimagnets and thereby three phases are found: cone, coplanar fan and a phase-separated one. As an 
application, we map out the phase diagram of a 3D helimagnet which consists of frustrated J1-J2 chains as a 
function of frustration and an interchain coupling. Moreover, we also calculate the stability of the 2-magnon 
bound state to investigate the possibility of the bound-magnon BEC. 
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I. INTRODUCTION 

Magnetic frustration introduces several competing states 
which are energetically close to each other and thereby desta- 
bilizes simple ordered states. One way to compromise two 
or more competing orders is to assume a helical (spiral) spin 
structure!. In this paper, we discuss the high-field behavior of 
a spin- 1/2 Heisenberg model with generic interactions: 

h = Y,j,,s,-s,+hY,s]. (1) 

For the simplest case with one magnetic ion per unit cell, one 
can easily find the classical ground state by minimizing the 
Fourier transform of the exchange interactions: 
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where the summation is taken over all j-sites connected to the 
i-site by Jij. When e(q) takes its minima emin at q = ±Q, 
helical order with the wave number Q or — Q appears (±Q 
are not equivalent to each other). 

When the external magnetic field is perpendicular to the 
spiral plane, the spiral is smoothly deformed into the so-called 
cone state (FiglT) and this persists until all spins eventually get 
polarized at the saturation field He. When the system has an 
easy-plane anisotropy and the external field is applied in the 
spiral plane, on the other hand, the system undergoes a (first- 
order) metamagnetic transition into a coplanar /an phase^. 

One of the simplest models which exhibit, at least in 
the classical limit, the helical order is a three-dimensionally 
coupled Heisenberg chains with nearest-neighbor- (NN) Ji 
and next-nearest-neighbor (NNN) J2 coupling. Because of 
strong quantum fluctuation in one dimension, the spin-singlet 
ground state of a single decoupled S ~ 1/2 chain can be 
quite different^ from its classical counterpart. However, it 
is generally expected that interchain couplings may eventu- 
ally stabilize the classical helical order. In fact, many com- 
pounds which contain these ID-chains as subsystems and 
display magnetic long-range orders are known (see, for in- 
stance, TABLE I. in ReflSI). Despite this naive expectation, 
even relatively mild quantum fluctuations in three dimensions 



may destabilizes the classical ground state in some frustrated 
systems^. Therefore, it would be interesting to explore the 
possibility that quantum fluctuation replaces the classical cone 
state with other stable ones e.g. a coplanar fan. 

Another interesting feature peculiar to the quantum case is 
that for a region slightly below the saturation field, we can 
view the emergence of various kinds of (weak) magnetic or- 
der as Bose-Einstein condensation (BEC) of magnons which 
enables us to use the full machinery of many-body theories^. 
The concept of magnon BEC has been successfully applied 
to explain various experimental result o^''°'^' . By using dilute- 
Bose-gas approach, Batyev and Braginskiil^ discussed mag- 
netic structure near saturation from a general point of view 
and concluded that this is the case if a certain condition for 
the bosonic interactions is satisfied. 

Recently, helimagnetism attracts renewed interest in the 
context of multiferroicitjii^ and multiferroic behavior has been 
reported for various helimagnets. For example, a helimag- 
netic material LiCuV04 may be viewed as coupled quantum 
S" = 1/2 J1-J2 chains and, as is expected from the classical 
theories, exhibits helical spin order'- and ferroelectricityi^iii 
simultaneously under moderate magnetic field. When the field 
is very high, on the other hand, this compound shows modu- 
lated collinear order-i^, which contradicts with the aforemen- 
tioned classical prediction^, and this suggests that quantum 
fluctuation plays an important role. In these multiferroic ma- 
terials, external magnetic field provides us with a way of con- 
trolling polarizationii and it is crucial to understand magnetic 
structures in high magnetic field. However, except for one- 
dimensional cases, only few reliable results are known for 
quantum systems so far. Our aim in this paper is to determine 
the stable spin configurations of 3D spin- 1/2 helimagnets in a 
fully quantum-mechanical manner. 

The organization of the present paper is as follows. In 
Sec. ini we describe how magnon BEC technique is used to 
determine possible magnetic structures slightly below the sat- 
uration field. Our dilute Bose gas approach predicts that in 
general there are at least three types of quantum phases (cone, 
fan and an attraction dominant phase; see FIG[T]l in the high- 
field region. By mapping to an effective Lagrangian, we study 
the low-energy properties of the cone and the fan phases in 
Sec.|III]and show that, on top of the standard Goldstone mode, 
there is a yet another gapless mode in the fan phase which cor- 
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responds to translation. 

In general, when Q is commensurate with the underlying 
lattice, lattice symmetry allows several higher-order interac- 
tions which may pin the above gapless translational motion. In 
Sec. lIVI we discuss the effects of commensurability on these 
phases. 

As a concrete example, the phase diagram of a model of 
coupled S = \/2 J1-J2 chains {J1-J2-J3 model) is consid- 
ered in Sec. [V] In this model, we also study the stability of 
2-magnon bound state. If the bound-magnon BEC occurs, the 
transverse magnetic moment vanishes. As a result, we find 
that, on top of the above three phases, a spin nematic phase 
appears. 




FIG. 1: (Color online) Two spin structures considered here: 'cone' 
(left) and 'fan' (right). In the fan structure, the spins are lying in a 
single plane (coplanar). The Q-vector is pointing along the c-axis. 



II. GENERAL FORMALISM 

A. Mapping to dilute-Bose gas 

To apply the powerful Bose-gas technique, we first rewrite 
the original spin model in terms of bosons. When the external 
field H is larger than the saturation field He, spins are fully po- 
larized along H (downward, here) and any spin-flip excitations 
from this reference state can be expressed exactly in terms of 
hardcore bosons. Specifically, we write the spin operators as: 



(3) 



Then, the original spin Hamiltonian in general may be rewrit- 
ten in the following way: 

(4) 



q,k,k' 



'^(q) = e(q) - en 

He = e(0) - Emin 



, M = He-H, 
= 2(e(q) + U) 



(5) 



where N is the number of lattice sites. In helimagnets, the 
single-spin flip excitation e(q) takes its minima Cmin at two 
inequivalent q-points ±Q. The external field H controls the 
chemical potential 1^1 and the on-site interaction 00) has 
been added to impose the hardcore constraint. 

In what follows, we consider a cubic lattice and assume 
that helical- and ferromagnetic/antiferromagnetic order oc- 
cur along the c-axis and in the ab plane, respectively (i.e. 



Q = (0, 0, Q) or (tt, it, Q)). Also, in order to avoid confu- 
sion, we use the indices (a, b, c) for the real-space coordinate 
and reserve {x,y, z) for the spin directions. 

We see that magnon BEC occurs when the external field is 
smaller than the saturation field: H < He (/J, > 0). Although 
the hard-core formulation is valid only for spin- 1/2, it can be 
generalized^^, with a little modification, to arbitrary spin-S*. 



B. Ginzburg-Landau analysis 

The thermal potential per site E/N of the dilute Bose gas 
is determined by the interaction among the condensed bosons 
at q = ±Q and the ground-state Boson densities p±Q are ob- 
tained by minimizing E/N. If we denote the renormalized 
interactions between the same bosons and that between differ- 
ent ones respectively as Fi and the energy density E/N 
is given by 



E 1,2 

N = 2^' (^Q 



P-q) + ^2 PQP-Q - P{PQ + P-q), 



= \iTi+ T2){PQ + p-Qf + T^Kpq - P-q) 



- p-Ipq + p-q) 



(6) 



where pq = \{Pq)\'^/N. First we note that the energy func- 
tion E/N has discrete Z2-symmetry Q <-> — Q. When the 
external field H is sufficiently close to the saturation field He, 
we may expect that the condensed boson is dilute and one can 
safely use the ladder approximation-^*' to calculate the inter- 
action vertex (see Fig|2]i: 



r,(k,,k2) 



Fq,(kl,k2)Vq_c 



N ^cj(ki+q') + cj(k2-q') - t^(ki) - tj(k2) 



(7) 



From this, one obtains the parameters Fi and F2 as^ Fi 

ro(Q, Q) , r2 = ro(Q, -Q) + r_2Q(Q, -Q). 



ki+q k, 



rq{k,.kj)= iq + 



k2-q ^2 



ki+q' ki+q 



r„(ki,kj) 



iq-q' 



krQ k2-q 



FIG. 2: Ladder approximation to the interaction vertex. 



1. Cone phase 

Different phases appear according to the values of Fi 2. If 
F2 > Fi > 0, the ground state is given by pq = p = /i/Fi, 
p_Q ~ (or vice versa) and E/N = —fi^/ (2ri). Hence, the 
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spin configuration is determined as: 

(A) - Vpexp{±*(Q-Ri+0)} , (5f ) = -\+P^ 

(5f) = V^cos(Q-Rz + 0), (5f)=TVPsin(Q-R, + 0). 

That is, the cone state (the left panel of Fig[T]l, which exists 
akeady in the classical case?, is favored for r2 > Fi. 

It is easy to see that this phase exhibits the multiferroic be- 
havior. According to the so-called spin-current mechanism, 
a microscopic electric polarization P^j, which is associated 
with a pair of sites i and j, is given b y'^i''* 

P,, =7?e,, X ((S,) X (S,)) , (9) 

where 77 is a constant. When the external field H is paral- 
lel to the a(or &)-axis, the Q-vector is in the spiral plane, 
which is perpendicular to H, and the summation of the lo- 
cal polarization P, j over the lattice yields a finite polarization 
(X pQ sin Q parallel to Q x H. When H is along the c-axis, 
on the other hand, the local polarization sums up to zero and 
the system shows no ferroelectricity. 

2. Fail phase 

If Ti > and T = Ti + > 0, on the other hand, the 
two modes condense simultaneously and the ground state is 
determined as: pq = p^Q = p' = fi/T, ^ = —p?/T 

(A) = Vp' je'^'^-^'+'^i) _,_gj(-Q-R,+e2)| ^ 
(5f> = -^ + 4p'cos2(Q.R, + ^l^), (10) 

{S^) = 2V7cos(Q-Ri + A^)eT'^ . 

The two parameters Oi and 62 characterize arbitrary phases of 
the two condensates (/3q) and (/3_q), respectively and lead 
to two different low-energy excitations. Since {Sf) / {Sf) = 
— tan ^ii^' the spins assume a coplanar configuration (fan) 
shown in the right panel of FIG[T] 

The ferroelectric property of this phase can be seen again 
from eq.®. If one moves from one site to the next along the 
c-axis, spins change their direction periodically within a basal 
plane specified by the azimuthal angle {61+62) /2 (see FIGlU). 
Although the vector chirality on each bond ((S;) x (S^+eJ) is 
always pointing a fixed direction perpendicular to the basal 
plane, it changes the sign within a period; for the first half pe- 
riod, it is positive and for the latter negative. Hence the local 
polarizations Gc x ((S^) x (S^+e^)), when summed up along 
the c-axis, exactly cancel out and yield zero macroscopic po- 
larization (note that only bonds parallel to the c-axis give non- 
zero contribution). 

Here we would like to stress that the fan state here does 
not require any kind of easy-plane anisotropy and should be 
distinguished from its classical counterpart which exists only 
in easy-plane helimagnets^. It is interesting to observe that in 
the second case (fan) the ordinary superfluid order ((S*^) ^ 0) 
and the spin-density wave, where Sf modulates with momen- 
tum 2Q, coexist. 



3. Attraction-dominant phase 

When Fi < or Fi + r2 < 0, low-energy bosons around 
q = ±Q attract each other. If the energy (|6| is taken literally, 
first order transitions may be expected on general grounds. In 
some cases, this scenario may be the case and, on physical 
ground, we may expect bosons to "collapse" in real space. It 
might well be that as a subsequent phase a cone or fan phase 
appears via first-order transition. However, eq.® is based on 
the assumption that magnon BEC occurs in the single-particle 
channel and may not work when we expect magnon bound 
states stabilized by strong attraction. In fact, this conditions 
for Fi 2 implies nothing but instability in the one-magnon 
condensates. In Sec|V] we calculate the energy of the 2- 
magnon bound state in the concrete model. As a result, we see 
that the bound state tends to be favored in attraction-dominant 
phase. 

We summarize the (Fi, r2))-phase diagram in FIG. [3] In 
section [V] we shall calculate Fi 2 for a specific model and 
show that all three possible phases appear. 




FIG. 3: (Color online) Phase diagram in (Fi, r2)-plane. When 
Fi < or Fi + r2 < 0, the energy function E/N is, at least within 
single-particle BEC, unstable (unstable (i) and unstable (ii)). The 
phase 'unstable' (i) and (ii) are respectively characterized by formal 
solutions (p±q ^ 00, p=pq = 0) and p±Q 00. 



III. LOW-ENERGY PROPERTIES 

A. Effective Lagrangian 

Thus far we have described the general resultsi^. Now 
we study the low-energy properties of the two phases more 
closely. To this end, it is convenient to introduce two indepen- 
dent low-energy modes a and b through 

A -e''5-^'a(R) +e-''Q-^'6(R) . (11) 

The mass parameters corresponding to these modes are ob- 
tained from the low-energy dispersion. Defining k = q— Q, 
we may expand e(q) = e^nm + hkj / {2mij) + • • •, = 
a, b, c) where the summation over repeated indices is implied. 
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The symmetric mass matrix my can be diagonalized to give 
a standard form of the dispersion kikj/2mij = k'^ / {2m'-) = 
eg(k'). We assume > (for all i) for the stability of the 
minima at q = ±Q. Below we omit the prime over rrij for 
simplicity. As a result, we write down the following effective 
Lagrangian with the renormalized interactions: 



Leff — 



dh 



-{a*dta-adta*)- 



\dja\' 
2m,- 



+lia*a 



2m i 



^i\af + \bf)^T,\am' 



We note that this effective Lagrangian has U(l)xU(l) sym- 
metry; one comes from the axial (around the external field) 
symmetry and the other from an emergent translational sym- 
metry which does not exist at the level of the lattice. When 
Ti = this U(l)xU(l)-symmetry gets enlai-ged to U(2). 
On the fully saturated ground state, both a and b bosons have 
the energy gap —fi. 



where we have introduced 9u = (6*1 + 92)1^/2, 9v = (^i — 
92)1^/2. Now the meanings of the two angular variables ap- 
pearing in eg. ([Toll are clear; the field 0„ corresponds to the 
Goldstone mode associated with the spontaneous breaking of 
rotational symmetry in the x-y plane, while the other 9„ de- 
scribes the translational motion of the fan along Q (phason of 
the spin (S^-) density wave). The excitation spectrum of the 
^„-mode is readily obtained from ( fTSl ) as 



(12) a.(k) = ^2(ri +r2)p'eg(k) = ^2/xeg(k) aone(k)) . 

(16) 

This gapless excitation does not exist in the fan phase appear- 
ing in classical models with easy plane anisotropy^. Similarly, 
the excitation spectrum related to 9^, is given by 



a„(k) = y'2(ri -r2)p'eg(k) = ^2^;i-^Aieg(k). (17) 

At the transition point Fi = r2 from the cone to the fan, the 
phonon velocity of the 6*1, -mode vanishes indicating an insta- 
bility in the translational mode. 



B. Cone phase 

In the case of r2 > Fi > 0, the cone phase appears. Only 
one of the low-energy bosons (say, a) condenses and we pa- 
rameterize it as: a = + Sp e*^ with Pq = p = P-/^i- If 
we integrate out the massive Sp fields in the path integral, we 
obtain the following effective Lagrangian: 



-(b*dtb-bdtb*) 



(Ml! 

2Fi 



2m. 



2m 1 



+ ^{l--±)b*b-^\b\' 



(13) 



From this, we can read off the excitation spectrum of the 0- 
mode as: 



f^cone(k) = y2^eg(k) , 
and that of b field acquires a gap — 1). 



(14) 



IV. EFFECTS OF COMMENSURABILITY 

In this section, we consider effects of commensurability on 
the ground state. Since our system is defined on a lattice, any 
types of interactions which are allowed by the symmetry may 
be added to the effective Lagrangian ( fT2] i. Specifically, we 
require invariance under 



(i) global U(l): (a, 6) 

(ii) lattice translation: 
{5: lattice period) . 



\a,b) 



b^b^-'"^-^ (18) 



For generic incommensurate values of Q, the effective La- 
grangian ( fTSl ) correctly describes the low-energy physics. If 
Q is rational (i.e. Q = Qo ^ irn/l: l,n are coprime), on the 
other hand, the following terms in general appear to break the 
translational U(l) symmetry explicitly: 



C. Fan phase 

The low-energy spectrum of the fan phase (Fi > F2, 
Fi + F2 > 0) exhibits quite a different behavior, since 
two bosons a and b condense simultaneously and the above- 
mentioned U(l)xU(l)-symmetry plays a crucial role. To see 
this more clearly, let us parameterize a = y/p' + Spi e**\ 
b = \/ p' + &P2 e**^ and integrate out the massive ^pi;2 fields 
in the path integral. Then Lgff reduces to: 



fan 



+ 



l2(Fi+F2) 
{dt9,? P 



2(Fi - F2) 2m^ 



2m 



■{dj9,f 



(15) 



-F3(at'6' + 6t'a') ~ F3/ cos /(6'i - 6I5 
= F3p' 003(72/ 6l,„) 



(19) 



If we treat the problem in a classical manner, we can imagine 
an infinite sequence of crystalline phases with superfluid order 
(spin analogue of supersolids). 

However, if we take into account quantum fluctuation, this 
devil's staircase structure is destroyed and only a finite num- 
ber of commensurate phases survive^!. To see this explicitly, 
we redefine the boson operator as: 



a(R) y-^ e-*'^<^-^a(R) , 6(R) ^ e+**'<^-^&(R) 
(JQ = (0,0,Q-Qo)) . 



(20) 
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Then, the effective Lagrangian is given as, 

\{-idj-5Qj)a\'' 



-{a*dta — adta*) 



2m i 



+ fia a+ - [b Otb — botb ) ^ fib b 



2m, 

^i\a\^ + \b\^) - r2|a|2|6|2 -^^(a*'6' + 6*'a0 



(21) 



From now, we concentrate on the fan phase where both a 
and b condense. As before, we integrate out the massive (5pi,2 
fields in the path integral, and ignore the terms which do not 
matter for low-energy physics when the superfluid density p is 
dilute and I > 3. For I = 2 case, term is the same order as 
Fi and F2 term in E/N and therefore fan phase appears for— 
Ti + iTal > F2. Thus, our approximation may not be justified 
for Z = 2 case. Now, Lquo reads. 



em 



l2(Fi+r2) 



2m 



2(Fi 



2m, ' ' 



{dj9y - V26Qj)^- Fap' cos V2W^ 



(22) 



Classically, if SQ is small, the Fa-term seems to pin the trans- 
lation mode 9y at the expense of the elastic energy. If 9y is 
pinned, on the other hand, the gapped zero-point fluctuations 
around the pinned value yield the (positive) quantum correc- 
tion to the ground state energy and thereby a soliton lattice 
with gapless excitations^ may be favored. From eq.(14) in 
Refl2U an incommensurate soliton lattice is stable for any SQ 
if the following inequahty is satisfied: 



r > Min 
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mi{Ti - F2) 



1/2- 



(23) 



where Min[- • • ] means that the minimum value with respect 
to i = a, 6, c should be taken. Therefore, at least in the dilute 
gas limit, i.e., just below the saturation field, a commensura- 
bility locking does not occur for I > 3. When p grows further, 
eq.(l23Tl may be violated for some small commensurability I 
and the locking occurs; the pitch Q is locked to its commen- 
surate value Qo until SQ exceeds the critical value^' 



6Q^, = (8A^)2TO,|F3|p'-Hl-aVl6)V"^c(ri - F2)M^ 

(24) 

It is interesting to see that the first term, which is obtained by 
classical calculation, has the same character as a classical fan 
phase in an easy plane, which has a width proportional to2i 
5Q^ cx |H-Hc|'-^ 

Before concluding this section, we would like to give a re- 
mark on the validity of our treatment. Above discussions as- 
sume the dilute-gas limit, where the scattering length is much 
smaller than the average interatomic distance p~^l^ . Specifi- 
cally, our approximation is valid when V i{mam]^mcpY^^ ^ 
1 is satisfied for j = 1 or 2. 



V. COUPLED J1-J2 MODEL 
A. Phases of a single Ji - J2 chain 

Before presenting our results for a 3D model (J1-J2-./3 
model), let us briefly review the known results for the 5 = 1/2 
J1-J2 chain (the case with J3 = 0) and discuss the connection 
to the phases found in Sec|II] In the case Ji > 0, near satura- 
tion, two dominant phases are found-^: (i) 'chiral phase (VC)' 
with finite vector chirality parallel to the magnetic fielda^ia^ 
and (ii) 'TL2' phase where the system is described by two 
Tomonaga-Luttinger (TL) liquids*^*^. Obviously, the former 
turns, after switching on an interchain coupling, into the cone 
phase. A close inspection of the two gapless TL modes near 
saturation tells us that the TL2 phase should evolve into the 
fan phase in three dimensions where we have two Goldstone 
modes. Yet another dominant phase 'TLl', for which a single- 
component TL gives a good descriptio n"^'^^ , is located in a re- 
gion where we expect a more conventional single-component 
BEC at Q = (0, 0, tt) ( J3 < 0) or (tt, tt, tt) ( J3 > 0). 

The ferromagnetic side Ji < is much more subtle as 
we expect BECs of rt-bound magnon states (n > 2) to oc- 
cur. In one-dimension (J3 = 0), on top of the VC phase 
described above, various phases related to bound n-magnons 
(2 < n < 4) have been found^^; (a) TL phases of 2- 
magnon bound states 'nematic' and 'SDW2', whose domi- 
nant correlation occur in respectively superfluid- and SDW 
channel, (b) 3-magnon TL 'triatic' and 'SDW3' (the mean- 
ings of them are evident) and (c) 'quartic' corresponding to 
4-magnon bound states. Our dilute-gas analysis predicts that 
inside the domes of the attraction-dominant phase (phase (iii) 
in FiglU one-magnon BEC becomes unstable toward various 
kinds of magnon bindings as has been discussed in Sec lIlB 3l 
and we may expect that the above n-magnon-based phases 
correspond to the attraction-dominant phase. We also study 
the stability of the 2-magnon bound state by the traditional 
approach2^ and discuss in later part of this section. 



B. 3D phase diagram 

Having established the formalism, we now consider a frus- 
trated spin- 1/2 model on a simple cubic lattice whose Hamil- 
tonian is given by 

-f^ = ^^{^lSr-Sr+g_, + J2Sr-Sr+2e<; +'^3Sr-Sr-(_ei} i (25) 
v,i—a,b 

where we label the three crystal axes by (a, b, c) and the spiral 
vector Q is pointing the c-direction. The J1-J2 chains are 
running in the c-direction and J3 controls the coupUng among 
adjacent chains. 

If we replace the spin-l/2s by hardcore bosons, we obtain 
the bosonic Hamiltonian (|4| with e(q) given by: 

€{q) = Ji cos qc + J2 cos 2qc + J3(cos qa + cos gh) . (26) 

The mass parameters are given by TOq — mi, = 1/| J3I, and 



1 / (4^/2 — 0^1 /4 J2 ) . The wave number Q characterizing 
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the condensate is given either by Q = (0, 0, Q) ( J3 < 0) or 
by Q = (7'',7r,Q) (J3 > 0) where Q = arccos(— J1/4J2). 
We solved eq.(|7]i to determine the spin structure of our J1-J2- 
J3 model (see Appendix lAl for the details). As a result, we 
obtained the phase diagram shown in FIG|4] 
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FIG. 4: (Color online) Phase diagram slightly below saturation 
(H < He) mapped out in ( Ji , J3) -plane obtained from the one- 
magnon-BEC approach (J2(> 0) is used to set the energy unit). 
Note that only the region —4 < J1/J2 < 4, where cone structure 
with incommensurate Q is expected classically, is shown. Two Bose- 
condensed phases (i) cone phase and (ii) coplanar fan phase as well 
as the phase-(iii) which is characterized by Fi < or Fi + r2 < 
are shown. The region | J3I/J1 ^ 1 is omitted (see the text). Inset: 
The same phase diagram for the large negative interchain coupling 
(J3 < -0.1). 



On top of the ordinary cone phase, strong quantum fluc- 
tuation in = 1/2 systems stabilizes two new phases: the 
coplanar fan ((ii)) and the attraction-dominant phase ((iii)). In 
the phase-(ii), both gauge symmetry and translation symme- 
try are broken simultaneously and as a consequence we have 
two different low-energy (Goldstone) modes r2 u(k) (eq.dTSIl) 
and r2u(k) (eq.dTTl)). In the phase-(iii), strong attraction may 
imply instabilities toward other phases e.g. conventional fer- 
romagnetic one or more exotic multipolar ones^. For J3 0, 
low-energy quantum fluctuation destabilizes F and our ap- 
proach cannot be extended to J3 = continuously (F be- 

1 /2 

comes 0( J3 ) and Fi F2 at the leading order in J3). 

For the ferromagnetic J3, the cone phase (region-(i)) gets 
wider and wider and the boundary between the cone- and the 
phase-separated phases approaches the classical phase bound- 
ary J1/J2 = —4 as IJ3I is increased (see the inset of FlGlH). 
This is easily understood since for very large negative J3 all 
spins sitting on each a6-plane behave like a single large spin 
to which classical analysis is applicable and the system may 
be thought of as a single chain running in the c-direction. For 
the antiferromagnetic coupling (J3 > 0), these novel phases 
((ii) and (iii)) appear only in the weak-coupling ( J3 <C J1.2) 
region. 

To see the possible magnon binding more clearly, we plot 
Fi in Fig|5] Although Fs behave regularly in the most part 
of the phase diagram, Fi has poles on the boundary between 
the fan phase (ii) and the phase-(iii) as is seen in Fig|5] This 



implies that near the boundary between the phase-(ii) and (iii) 
the interaction among bosons becomes singularly large which 
may lead to new phases. Actually, a pole of a interaction be- 
tween two particles in general imply an existence of stable 
bound states. Thus, the one-magnon-BEC approach is not suf- 
ficient to see the ground state near this boundary. 
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FIG. 5: (Color online) Fi and F2 for J3/J2 = -0.01 plotted as a 
function of J1/J2. Also plotted are Fi of classical ID-chain (J3 = 
0), where we use the spin-wave expansion around the saturated phase 
and retain only the leading term of the 4-point interaction in S as in 
Ref0. In the classical case, the cone phase is always stable except 
at J1/J2 = 2. In the quantum case (J3/J2 = —0.01), Fi has two 
poles at the boundaries between the phase (ii) and (iii) in FIG|4l 



To highlight this point, we study the instability of the 2- 
magnon bound state. On the fully saturated ground state 
I FA/), the wave function of the 2-magnon bound state is 
given by j i^ii, j)S^ S^\FM), and the energy of this 
wave function can be exactly obtained by solving the two- 
body Schrodinger equation^"^. If the gap of the bound state 
closes earlier than that of the one magnon, the bound-magnon 
BEC will occur, and the nematic order emerges in the trans- 
verse direction. We show the region of the stable bound state 
in Fig|6l As a result, the nematic phase completely masks 
the fan phase which would have appeared on the ferromag- 
netic side Ji < 0. At a rough estimate, for —2.7 < J1/J2, 
the lowest-mode of the bound state is commensurate and for 
J1/J2 ^ —2.7 the one is incommensurate in the same way as 
in the ID J1-J2 chain^. 

Detailed results on the bound magnons will be reported 
elsewhere^. 



VI. SUMMARY 

By using the dilute-Bose-gas technique, we studied the 
high-field magnetic structures and low-energy excitations of 
three-dimensional quantum (S = 1/2) helimagnet. The 
method used is asymptotically exact when magnetization is 
close to saturation (we gave a criterion of 'proximity' in the 
end of section IV) and enables us to obtain reliable results for 
three-dimensional frustrated systems. 

In Sec.|II] we discussed various phases emerging from the 
BEC of hard-core magnons slightly below the saturation field. 
Although only the cone phase is expected in the classical 
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FIG. 6: (Color online) The same phase diagram as in Fig|4]when 
the 2-magnon bound state is taken into account. The phase (iv) is 
characterized by the condensation of the 2-magnon bound state and 
leads to the nematic order in the transverse direction. 



heHmagnets^, quantum fluctuation can stabilize the fan or at- 
traction dominant phases as well according to the renormal- 
ized interactions (Fi, r2). 

Then, the low-energy excitations of the cone and the fan 
phases were studied in Sec. Hill The hallmark of quantum he- 
limagnets is that one has two low-energy modes at k = ±Q 
and the low-energy physics is described by the effective La- 
grangian with U(l)xU(l) symmetry; one comes from the ax- 
ial (around the external field) symmetry and the other from an 
emergent translational symmetry. In the cone phase, only one 
of the two bosons condenses and there is one gapless Gold- 
stone mode. Meanwhile, the fan phase breaks both symme- 
tries and has two types of gapless Goldstone modes. 

In Sec.|IV] we considered effects of commensurability on 
the helical modulation vector Q. If Q is rational (i.e. Q = 
Qo = T^n/l: l,n are coprime), additional interactions appear 
in the effective Lagrangian and favors the gapped commen- 
surate phase. By examining the correlation due to the gapped 
zeropoint fluctuations, we found, for I > 3, that slightly below 
the saturation field quantum fluctuation destroys the commen- 
surate order 

We studied a concrete three-dimensional 5* = 1/2 model 
(J1-J2-J3 model) in Sec. [V] and mapped out its (high-field) 
phase diagram in Figs|4] and |6l An analysis assuming the 
single-magnon BEC predicts the existence of a fairly large 
region where single-magnon BECs may be unstable (phase- 
(iii)) as well as the cone and the fan shown in FiglT] To get an 
insight into the nature of this 'phase', we consider the possi- 
bility of a BEC of two-magnon bound states, which leads to 
the nematic order in the transverse direction. In fact, in a large 
portion of the phase-(iii) (a region marked as (iv) in Fig|6]l, we 
have a stable two-magnon bound states which condense first 
when the external field is decreased. In the ID J1-J2 model 
(i.e. J3 = 0), it is known that one has multiple-magnon bound 
states up to four-body^i and some parts of the new phase {ne- 
matic (iv) in Fig|6]) could be replaced by the condensed phases 
of these bound states. 



Finally, we comment on the relevance of our study to real 
systems. Interests in multiferroicity sparked an intensive 
study of various helimagnetic materials, among which one 
can find many examples of coupled J1-J2 chains. For ex- 
ample, LiCuV04 is characterized by edge-sharing CUO2 pla- 
quettes and may be modeled by the S* = 1/2 J1-J2 chain 
with negative Ji. Neutron diffraction and ab initio calcula- 
tions suggestedi^, as well as Ji and J 2, various kinds of in- 
terchain interactions J^, ..Jq. Although the stacking of J1-J2 
chains is different from what is assumed here, our method can 
be readily generalized to include more realistic cases and we 
hope our approach will shed some light on magnetism of these 
quantum helimagnets. 
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APPENDIX A: HOW TO TREAT THE LADDER DIAGRAM 

For convenience, we briefly summarize the method of cal- 
culating r. To obtain Fi = rq=o(ki = Q,k2 = Q) and 
r2 = ro(Q, -Q) + r_2Q(Q, -Q), we solve the following 
integral equation in the case of ki 2 = ±Q (ci'(ki 2) = 0): 

r (U \c \-V ^ \^ rq.(kl,k2)Vq,q^ 

(Al) 

where = 2(e(q) + U). In what follows, we do not write 
the argument ki_2 of F explicitly, and denote ( ). 

Since (e) = 0, we sum up the both side of ( lAlb with respect 
to q and obtain. 



(F) = 2J7 I 1 - — V — 

^ ' 1 N ^ cj(ki +q' 



q' 



Using this equation, jAlb is simphfied to, 

Fq = 2e(q) + (F) 

M 2^ 



(A2) 



(A3) 



A^^a;(ki+q')+w(k2-q') 



Additionally, if we assume the limit U 00, eq.( IA2l i reads: 



1 V- 



q' 



N ^ w(ki + q') + tj(k2 - q') 



= . (A4) 
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Now, the problem is reduced to solve ( IA3I ) and ( IA4I ) simul- 
taneously, which are free from the infinite term U. Next, we 
expand Fq in lattice harmonics. Since 



We note that (F) and Ai are independent of q, but depend on 
(ki, implicitly. If we substitute this, eq. dASt reduces to. 



e(q ~ q') —Ji (cos Qc COS q'^ + sin qc sin q'^) 

+J2 (cos 2qc cos 2q^ + sin 2qc sin 2q^) 
+ J3 (cos cos q'a + sin ga sin g', 
+ cos 5f, cos q'l, + sin sin g^) , 



(A5) 



we introduce 



Fq = (F) + JiAi cos qc + J1A2 sin qc + J2A3 cos 2qc 



J2A4 sin 2(jc + J3A5 (cos ga + cos qi,) 



(A6) 



2 Ji cos g' ^ „ T \ / , , 1 

-—T^,-2J,\cosqc+\jiA2 + -J2 



N ^ tj(ki + q') + cj(k2 - q') 



2 Ji sin g^ 



iV ^ Lj(ki + q') + w(k2 - q') 



-Fq/ singe 



2 J2 cos 2g;^ 



iV ^ tj(ki + q') + cj(k2 - q') 



-Fq, - 2J2 cos2gc + J2A4 



■E 



2 J2 sin 2g^ 



TV ^ u;(ki + q') + tj(k2 - q') 



— Fq/ sin2gc 



(-. , 1 ■^-^ >/3(cosg' + cosgH , , 

-^3^5 + ^ E \ ,\ \ ^A ^q' - cosg, + cosgfc = , 

iV ^ w(ki + q') + cj(k2 - q'j ^ / 



(A7) 



where we use the relation 



-E' 



sm q. 



N ^ w(ki + q') + cj(k2 - q' 



TrFq/^O 



(A8) 



To satisfy the eq. dAVI ) for arbitrary q, the coefficients of 
trigonometric function of q must be 0. For convenience, we 
define 



nj{ki,k2) ^^Yl 



Uq')Uq') 



N ^ uj{ki + q') + Uj{k2 ~ q') ' 



where 



T(g) = (1, cosgc, singe, cos2gc, sin2gc, cosga + cosgb) . 

(AlO) 

(A9) Then, (|A4] | and dATl i are put together into 





Til 


J1T12 


J1T13 


J2T14: 


J2T15 


J3T16 




/(r)\ 






2t21 


1 + 2 Jir22 


2J1T23 


2^2^24 


2J2T'25 


'2J3T26 


Ai 




2 




2t3i 


2J1T32 


1 + 2 J1T33 


2 J2T34 


2 J2T'35 


2J3T36 




A2 









2t41 


2J1T42 


2J1T43 


1 + 2 J2T44 


2 J27'4.5 


2J3T46 




A3 




2 




2t5i 


2J1T52 


2J1T53 


2J2T54 


1 + 2 J2T55 


2J3'''56 




Ai 

\aJ 







I 


T61 


J1T62 


JlT63 


^2T64 


J2Te5 


1 + J1T66 


) 







(All) 



This equation can be solved by calculating Ty numerically. If symmetry and we obtain, 
we evaluate at (ki, k2) = (Q, Q), A2 = A4 = due to the 



Fi = (F) +J1A1+ J2A3 + 2J3A5 



(A12) 
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If we evaluate at (ki , k2 ) = (Q, —Q), we obtain The following procedure is the same as in the former case, and 

r2 is given by, 

r2 = [(r) + JiAi + J2A3 + 2J3A5] 

+ [(r) + JiAi cos(-2g) + J1A2 siii(-2Q) 
+ J2A3 cos(-4Q) + J2A4, sin(-4Q) + 273^5] , 

(A13) 



Although above we review the straightforward method, we 
can calculate r2 more simply if we introduce, 

rq-Q(Q,-Q) + r_q_Q(Q,-Q) 

=2(r) + Ji cos + J2 ^2 COS 2qc (A14) 
+ JzA''i{cosqa + cosqb) . 



T2 = 2(r) + JiA;cosQ + J2^2Cos2Q-2|J3|Ajj . (A15) 
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